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Abstract 

The heterogeneity of composite leads to the extra charge concentration at the boundaries of 
different phases that results essentially nonzero effective electric susceptibility. The relation be- 
tween tensors of effective electric susceptibility Xef an d effective conductivity a e j of the infinite 
two-dimensional two-component regular composite with rhombic cells structure has been estab- 
lished. The degrees of electric field singularity at corner points of cells are found by constructing 
the integral equation for the effective conductivity problem. The limits of weak and strong contrast 
of partial conductivities are considered. The results are valid for thin films and cylindrical 

samples. 

PACS numbers: 73.25. +i, 73.40.-c, 73.40. Jn, 73.50.-h, 73.61.-r 
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I. INTRODUCTION 



The evaluation of effective properties for two-dimensional (2D) two-component com- 
posites, which determine the behavior of the medium at large scales, given rise by Kelleilll 
and Dykhndl, remains a topic of high activity. Among different approaches (variational 
boundsi@; asymptotici'!; numerical^; network analogue!'!) used to consider this problem, 
the analytical approach, being classical problem of mathematical physics, is surprisingly 
very difficult. Exact values of effective parameters are of great interest even though these 
values are established in idealized models. It seems that explicit formulae are available 
only as exceptions. Such formulae which solve the field equations were obtained for two- 
component regular checkerboard with squared, rectang ularSS and triang ulaiH unit cell 
using complex-variables analysis. Another technique (integral equations) was used in the 
recent papers dealt with squared and triang ulaJl! regular checkerboard. 

Almost all these studies were directed towards effective conductivity a e f evaluation de- 
spite of important fact that the heterogeneity contributed to the conducting composite some 
dielectric properties. The homogeneous metal does not possess static dielectric properties 
(such as electric susceptibility) because only core electrons can contribute there, but their 
influence is obviously small. However, heterogeneity leads to the extra charge concentration 
at the boundaries of different phases, which results essentially nonzero effective electric sus- 
ceptibility Xef- The implication follows that the relation between the conductivity a e f and 
susceptibility Xef effective tensors must exist. 

In the present paper we will consider the regular 2D two-component rhombic checker- 
board and derive such relation. This middle-symmetric structure belongs to p'jnm -plane 
grouplll and gives rise to anisotropy of a e f. In some sense this anisotropic model is more uni- 
versal than the regular 2D two-component rectangular checkerboard (c'mm -plane group). 
Really, the effective electric properties are mostly determined by the corner points of the 
cell, where the electric field is singular!!. The structure of composite near these points in 
rhombic checkerboard is governed by arbitrary angular variable. 
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II. INTEGRAL EQUATION 



The regular checkerboard structure is composed of rhombic conducting cells with isotropic 
homogeneous conductivities <j\ and a 2 , hereafter <j\ > o 2 . The backbone of such structure 
can be represented as the set of images of the letter "X" with infinitely long legs, which 
are shifted up and down to distance 2Acos |r, N = 0, ±1, ±2, . . .; a is the smallest angle 
between legs (see Figure |1]). The side of the cell is scaled by unit length. Such arrangement 
of the checkerboard allows to generate the kernel for the integral equation by single series 
summation rather than by double summationBIl!. 

We will consider the external unit field E be applied in vertical direction Y. It is one of 
the principal axes of the effective material tensors. Another axis X may be considered just 
changing a — ► tt — a. 

Let us proceed with solution of Laplace equation for a scalar potential 0(r) at the infinite 
plane 5* 

0(r) = -E y - Air f G(r,r 1 )p(r 1 ) d 2 r x , G(r, rj.) = — In |r - r x | , (1) 
Js Aix 

where G(r, ri) is the two-dimensional Green function and p(r) is a charge distribution at 
the plane. The boundary conditions at the edge relate normal components of the field E n 
and of the current density j n 

- = Anp(t) , 3n (t) = „ X E® = a 2 E® , ( 2 ) 

where a new variable t is introduced to measure the distance along the edge of a unit cell 
counted from the cell corner and p(t) hereafter is the charge distribution at the edge. The 
boundary conditions (0) allow to write master equation 

E^+E^ = A -fp(r), Z=*^*, 0<Z<1. (3) 

Finding the corresponding derivatives E$ (see Appendix [S]) we come to the integral equation 

^l^p(t) = sin | - 4 , dt' , (4) 

where a new function g(t) reflects a periodic interchange of the constituents (<7i and cr 2 ) 
with variation of argument t 

0(t) = sgn[mod(t,2)-l], g(-t) = -g(t) ; g{t) = -1 for < t < 1 , (5) 
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the function mod(t, 2) gives the remainder on division of t by 2 and sgn[x] gives -1, or 
1 depending on whether x is negative, zero, or positive. The kernel Ki(t,t') is given by 
formula 

k k + t' 



k=— oo 



+ 



(t - t') 2 tan z + (t-f- 2k) 2 cot f (t + t') 2 tan § + (t - f - 2k) 2 cot § 
It is worth to represent the kernel for further summation as 

k k + t' 



(6) 



1 a + 
- tan — 
4 2 , ^ 

fe=— oo 



+ 



(k - kx)(k - k*) (k - k 2 )(k - k* 2 ) 
where the zeros k%, k\ and k 2 , k\ of both denominators in (|]) read 



k±, ki 



t-t' 



+ i 



\t-t'\ 



t-t' 

tan ^ , k 2} k* 2 = — - — ± i 



a 
2 ' 



\t + t'\ 



a 

tan- 



(7) 



(8) 



Making use of identity 



k + t' 



Im 



'h + t'\ 1 



(k — ki)(k — k*) "\k-kijlmki 
we can evaluate (0) in the sense of principal value and reduce essentially the kernel of integral 
equation @ 



7r f ki ■ cot 7rA;i (A; 2 + t') ■ cot 7r/c 2 

Al (t,t ) = lm < : : 1 : : 



(9) 



The further simplification of the kernel K\(t,t') can be continued by usage of trigonometry. 
Introducing p(t) = p(t)g(t) we obtain finally the integral equation 



2 , 1 a 

~ rp(t) = sin - + 

Z Ti 2 



+oo 



p(t')K 2 (t } t')g(t') dt' 



(10) 



where 



K 2 (t,t') 



tan s sin ir(t — t') — sinh ( 7r(t — t') tan ^ J tan ^ sin 7r(t — t') — sinh ( 7r(t + t') tan | 



cos7r(t — t') — cosh ^7r(t — V 



-f tan? 



+ 



cos 7x(t — t') — cosh (ii(t + t') tan ^ 



The function p(t) being a solution of integral equation ( pTJ]) makes it possible to find an exact 
expression of the effective conductivity tensor a e f (see Section fV\j . 



III. ASYMPTOTIC BEHAVIOR OF p(t) NEAR THE CORNERS 

We start this Section with two algebraic properties of the function p(t), which will be 
used in order to simplify further calculation. These are the parity and periodicity of the 
functions p(t), p(t), which are following from (|I0"D 

p(-t)=p(t), p(t + 2)=p(t) — p(-t) = -p(t), p(t + 2) = p(t). (11) 
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They are in full agreement with physics of the charge distribution pit) along the edges of 
the cells. A proof follows from an accurate evaluation of integral in fllPp. Indeed, the parity 
property follows due to fl5|) and identity K 2 (—t, —t') = —K 2 (t,t') 

9 1 rv r—oo r+oo 

-p(-t) + - sin £ = - / p(-t')K 2 (-t, -t')g(-t') dt ! = / p(-t>)K 2 (t, t')g{t') dt' . 

(12) 

The periodicity could be proven in the similar way. 

A similar integral equation was appeared in Ref. [14] for the two - component checkerboard 
with square unit cell. Its solution is presented by the means of Weierstrass elliptic function 
p sq (t) oc p K {t), where simiK = Z, and is found by inspection its behavior near the branch 
points t = and t = 1 

Eg>(t) ~ Psq (t) W ± , E®(t) „ Psq (t) ^ (1 - t)^ , A = Ax = k • (13) 

The equality of the exponents Ao = Ai is here essential. Already the two-component 
checkerboard with triangle unit celllll breaks the validity of (|T3|), that made unattainable an 
explicit0 solution, but only an efficient approximate method was proposed. 

The rhombic structure, discussed in the present paper, also gives rise to distinct expo- 
nents. Asymptotic behavior of p(t) near the branch points t — and t — 1 can be found 
from the equation ( [TD[ ) (see Appendix |A]) 

it) ~ p{t) ~° , sin 7cp = Z sin [a + ^ (tt - 2a)] , (14) 

E® (t) ~ pit) (1 - t) 2f11 , sin tt/xi = Z sin [a + p x (2a - tt)] . (15) 

Then 

1 71 7T 

l^o(Z) = (JLxiZ) — k — — arcsinZ if a = — ; ^o(Z) = —p,\{—Z) if a/- (16) 

tt 2 2 

and for a large contrast in conductivities a 2 <Ti, Z ~ 1: 2//o = 1, 2/ii = a/(7r — a). 

This shows that the generic rhombic cell (a^) does not lead to the solution of (|i~0|). 
which can be built out by simple rescaling of Weierstrass elliptic functions. An explicit 
solution of integral equation remains to be performed. 

It turns out that the integral equation (|T0D, obtained in the present Section, is sufficient 
to establish an exact relation between effective electric conductivity a e f and effective electric 
susceptibility Xef, which was not to our knowledge discussed earlier. 
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IV. EFFECTIVE SUSCEPTIBILITY OF RHOMBIC CHECKERBOARD 



Let us consider the polarization of the rhombic checkerboard at the scales large, compared 
to the size of the cells. The effective electric susceptibility is the tensor which defined by 
P = Xe/E where P is polarization and E is external electric field. In the reference frame 
(Figure |l|), when Xef is diagonalized, its y-component is determined by induced dipole 
moment d y per unit square: xtf = d y /S, where 

*v= E ^ = 2cos- E / (t + 2k)p(t)dt (17) 

over all Z fc=-00 °° 

charges qj 

is dipole moment of area S = L x L y and L x ,L y are the sizes of a sample. The summation 
covers all induced edge-charges qj which are placed within the area S. The sample which is 
composed of 2N X x 2N y unit cells has the area S = 2N X 2 sin | x 2N y 2 cos | = 8N x N y sin a. 
The formula ([[7]) can be reduced making use of parity and periodicity properties (|Il]) of 
p(t). Its accurate evaluation reads 

a r°° a n=Nx r 2n+1 a r 1 

d y = AN y cos — / tp(t)dt = 4N y cos- / tp(t)dt = 8N x N y cos- tp{t)dt . 

2 J— oo 2 n __jY J2n—1 2 J — 1 

Taking into account the inparity (|5]) of the function g(t) we obtain finally 

a /•! If 1 
d y = l6N x N y cos - / tp(t)dt , and x y ef = —— tp(t) dt . (18) 
2 Jo sm g Jo 

We define also the effective conductivity <r e / as a ratio of the current J = J j n (t) dt through 
the x-cross-section of the checkerboard per the unit length to the applied field E = 1 

< S = A \ 0l ° 2 f P{t) dt , ot f (a) = o* f {* - a) . (19) 
1 sm | (Ti — cr 2 Jo J J 

Due to KelleiS the principal values of the tensor a e j satisfy the duality relations 

a x ef (a) ■ a y ef (a) = a x a 2 . (20) 

Relating now two physical quantities ([18]), (|19D, we make use of auxiliary integral equation, 
obtained by integrating the equation (|T0|) (see Appendix 0) 

— / p(t)dt = —sm-- tp{t)dt. (21) 

(Ti — (T 2 JO 47T 2 JO 

The last relation could be rewritten in new notations (|TED, ( P^| ) 

AtixI f = 1 — > and similarly 4vrXe/- = 1 — . (22) 



One can think that a%, which appeared in p2|), breaks the universality of the formu- 
lae. Actually, the denominator contains the maximal value of partial conductivities 
a max = maxOi,cx 2 ). 

In fact, we have established the tensorial relation in any reference frame 

4vrx e / = / —&ef, (23) 

<^max 

where I is an identity matrix. Formula ( p3| ) results in the particular square-checkerboard 



case: = ^xtf = 1 - yW "!- 

Let us consider now two different cases of the weak and large contrast in partial conduc- 
tivities. 

1. (tri - <r 2 )/(Ti < 1, or Z < 1: 

In the first order on Z the integral equation ([H]) gives according to definition ( |T9| ) 

P® = ~ sin | — > = = a x (l - Z) (24) 

and therefore 

4^ = ^Xl f = Z . (25) 

2. cr 2 <C cti: 

In this limit the corner points t — > of the cell become important. Making use of the 
distribution (|T4| ) for the fields and charge ~ p(t) ~ ^—2^0 one can g n( j approximately in 



leading terms 



// = ~-^==W^ — ► ^r^V^' ( 26 ) 
/a(7r — a) V 1 



Actually, the coefficient was found inQ : A = ya/(n — a) cot (a/2). Formula ( Pq ) implies 
as well 



4^ f = l-A^, 4^ = 1-1 /^. (27) 
u o"i Ay G\ 



V. CONCLUSION 



1. We have derived the integral equation for the effective conductivity problem for the 
regular 2D two-component rhombic checkerboard. An asymptotic behavior of the electric 
field was investigated near the singular points t — and t — 1. 



2. The heterogeneity of composite leads to the extra charge concentration at the bound- 
aries of different phases that results essentially nonzero effective electric susceptibility. The 
exact relation ( p3|) between the two most important electrical properties, namely, effective 
conductivity a e f and effective susceptibility Xef, of rhombic composite was established. An 
absence of specific angular parameter a in this formula make us possible to conjecture its 
validity for any anisotropic two-component structure. It is shown that the tensor of electri- 
cal susceptibility has surprisingly simple structure in both cases of large and small contrast 
in partial conductivities a\, a^- 

3. The relation derived in present paper is definitely valid for cylindrical samples. It is 
also valid for thin films due to the conducting nature of the constituents, which confine the 
electric field inside the conductor. 
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APPENDIX A: DERIVATION OF INTEGRAL EQUATION ©. BEHAVIOR OF 
ITS SOLUTION NEAR THE BRANCH POINTS. 

We define the variables 

x = t sin — , y = t cos — , x' = ±t' sin — , y' = t' cos — \-2k cos — (Al) 
2 2 2 2 2 v; 



and normal vector to the edge 



n = ( cos — , — sin — ) . (A2) 



2 2. 

Here k is an ordinal number of the "X" image. Taking in mind the contributions from the 
both left (1) and right (r) edges of the rhombic tile we will find the derivative dcj)/dn = 
(n V) 

d(f) „ . a ^ roo l fdr 2 (t,t')\ ~ r°° //N I fdrf(t,t'Y 



s 



which lead after simple algebra to equation 

\ K } + 4 2) ) = -Eo sin | + 4 jH p(f )#i(f , ^ , (A3) 
where the kernel K\(t,t') reads 

„ / /n fc t -\- k 

= k hoo L(* - ^O 2 tan f + (t - t' - 2A;) 2 cot f (t + i') 2 tan § + (*-*'- 2A;) 2 cot f 

(A4) 

Taking now £?o = 1 we arrive at (||). 

Below we consider the asymptotic behavior of p(t) near the branch points t — and 
t = 1. 

• * — ► 0. 

Let us assume the power behavior p(t) oc | ^ | 2 a*o . S g n ^ anc j i 00 k f or the exponent p, . 
The main singularity comes from integral in ([10]) in the vicinity t' — > 0. The kernel K 2 
behaves as 



^(t-t') 2 + (t + t') 2 tan 2 f ' (A5) 



that gives the asymptotic behavior 



pit) 1 . f+°° dt' t' pit') 
- — sin a 



Z 7T 7-oo t 2 + it 1 ) 2 - It t'cosa 
Defining a new variable z = t'/t we obtain 



7T 



+ / ^^tttttt • (A6) 



Z sin a Vo 1 + 2 2 — 2z cos a Jo 1 + z 2 + 2z cos a 

The evaluation of the last expression is based on the primitive fraction expansion with 
further usage of standard integrals and gives finally ([14]) . 
• t — ► 1. 

Let us assume the power behavior p(t) oc |1 — t| 2Ml • sgn (1 — t) and look for the exponent 
p,\. It is convenient to define new variables t = 1 — t, r' = 1 + 1' and consider the vicinity 
of branch point r — > 0, so p(t) oc |r| 2Atl • sgn (r). In order to deal with a singular part of the 
integral equation ([10]) let us differentiate the last over t 

2 d pit) f+°° . n d K 2 (t,t') , d pit) n . , 2(/ , . 
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The main singularity comes from integral in (|A7|) in the vicinity t' — > 1, or r' 
the kernel behaves as 

dK 2 (t,t') r ,r'-»o 4 sin a t'(t + t' cos a) 

dt 7T (r 2 + (r') 2 + 2rr'cosa) 2 ' 

Defining a new variable i> = t'/t we obtain 

7r//i r +ca dv f 2 ^ 1+1 (l + vcosa) r+°° dv v 2 ^ 1+1 (l - v cos a) 
Jo (1 + f 2 + 2v cos a) 2 io (1 + v 2 — 2t> cos a) 2 



Zsina 

Evaluating the last integrals we arrive at (|H>1). 



0, where 



(A8) 



(A9) 



APPENDIX B: DERIVATION OF INTEGRAL EQUATION (M) 



Reminding that the equation (|T0| ) is written in the sense of principal integral value therein, 
we average this equation at the large interval [— M, M], taking afterward its limit M — > oo 



M r 2 . 1 a 
-pit) + - sm - 

Z 7T Z 



dt 



2M y_M 

where due to ( |Tl~l) the left h. s. reads 

2 r 1 



/•M /-4-Af 

/ dt p(t')K 2 (t,t')g(t') dt' , (Bl) 

J-M J-M 



2M J-M 



f 1 ~ 1 

/ p(t) dt -\ — sin 

JO 71 



(B2) 



while the right h. s. could be simplified. Indeed, let us represent the kernel K 2 (t,t') as 
follows 



K 2 {t,t') 



(tan f + cot f ) 


sin7r(t-f) (tanf + cotf" 


sin7r(t - if) 


cos ir{t — t') — cosh 


7r(t + t 7 ) tan f 


cos 7r(t — t') — cosh 


7r(t-t')tan§ 



+ 



1 a d i , 
— cot — — In < cos nit — t — cosh 
7T 2 1 



1 a d ( . .. 

— cot In < cos 7i it — t — cosh 

vr 2dt I V ; 



7i (t — t') tan 
71 (t + t') tan 



+ 



(B3) 



The last three terms do not contribute to integration of the kernel over t, while the first 
term implies 

I pM r-M 



2 f 1 _ 1 a 

- / p{t) dt + - sin - 

Z JO 71 Z 



2M J-M J-M 



p{t')K z {t,t')g{t') dt'dt 



(B4) 



where 
K 3 (t,t') 



sin Tit 



2 f 

7 r- = Re < cot — 

sin a cos vrt - cosh (ir(2t' + t) tan fj sin a I 2 



t + i(2t' + t)tan 



a: 



2 

(B5) 
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Continuing the integration of the kernel K^(t,t') we notice that the function 

r+M 

Fit') = lim / K 3 (t,t') dt (B6) 

M-*oo J-M 

is 1-periodic function: F(t' + 1) = F(t'), that follows from the structure of the kernel 
K^(t,t'). Evaluating the integral in ( |B6| ) we arrive at the following 

F(t') = 2(1 - 2t') for < t' < 1 . (B7) 

The final integral equation 

2 f 1 - 4) C P(t) dt + - sin - = 4 f 1 t p(t)dt (B8) 
V Z J Jo n 2 Jo 

leads already to ([2"T|). 
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FIG. 1: Regular rhombic two - component checkerboard under electric field Eq: unit cell (left) and 
basic variables for integral equation (right). The distribution of the charges is drawn in accordance 
with chosen inequality a% > o~i- 
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FIG. 2: The branch points exponents fJ-o(Z) (dashed line) and ^i(Z) (dot-dashed line) for rhombic 
unit cell with a = ir/3. The exponent k(Z) (plain line) for square unit cell is also presented. 
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